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Abstract 

New parameter-sets of the semi-realistic nucleon-nucleon interaction are developed, by modifying 
the M3Y interaction but maintaining the tensor channels and the longest-range central channels. 
The modification is made so as to reproduce microscopic results of neutron-matter energies, in 
addition to the measured binding energies of doubly magic nuclei including ^'^'^Sn and the even-odd 
mass differences of the Z = 50 and N = 82 nuclei in the self-consistent mean- field calculations. 
Separation energies of the proton- or neutron-magic nuclei are shown to be in fair agreement with 
the experimental data. With the new parameter-sets M3Y-P6 and P7, the isotropic spin-saturated 
symmetric nuclear matter remains stable in the density range as wide as p < 6/>o, while keeping 
desirable results of the previous parameter-set on finite nuclei. Isotope shifts of the Pb nuclei and 
tensor-force effects on shell structure are discussed. 
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I. INTRODUCTION 



As exotic natures of unstable nuclei such as the new magic numbers and the neutron 
halos are disclosed by experiments, microscopic studies based on the nucleon-nucleon (NN) 
interaction become even more desired in nuclear structure physics. While the fully micro- 
scopic A^A^ (and NNN) interaction is still too complicated to cover large volume of nuclei in 
the periodic table despite significant progress [ItiS], the semi-realistic A^A^ interactions have 
been developed 0, [sj by modifying the Michigan 3-range Yukawa (M3Y) interaction 0,3, 
which was originated from Brueckner's G-matrix at the nuclear surface and expressed by 
the Yukawa functions. The modification has been made so that the saturation and the 
spin-orbit (is) splitting should be reproduced within the mean- field approximation (MFA). 
Owin g to the recently developed numerical methods j8|-ll|, self-consistent calculations in the 
MFA iMi and in the random-phase approximation (RPA) [l5| have been implemented 
using the semi-realistic interactions. Among the advantages of the M3Y-type semi-realistic 
interactions is that they contain realistic tensor channels as well as correct longest-range cen- 
tral channels originating from the one-pion exchange, which have been pointed out to play 
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The semi-realistic 



significant roles in Z- or A^-dependence of the shell structure 
interactions are suitable to investigate effects of these channels within the self-consistent 
MFA and RPA H. ^ 

were adjusted 
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While the parameter-sets of the M3Y-type interactions in Refs. |5| 
to the data on the nuclear structure, some of them have been applied to the nuclear reac 
tions ll9i and to the neutron stars 



20[ as well. In studying structure of the neutron stars, 
density-dependence of the symmetry energy is crucially important [2l|. It has been pointed 
out that the symmetry energy at low density is significant in nuclear reactions: e.g. the 



charge-exchange reactions [19 



22[. The symmetry 



and the multi-fragmentation processes 
energy at low density may also affect the so-called pygmy dipole resonance in neutron-rich 
nuclei However, the symmetry energy, particularly its density-dependence, was not 

sufficiently reliable in the previous parameter-sets in Refs. jH, 18|| as pointed out in Ref. 24 1. 
giving rise to instability of the symmetric nuclear matter at the density p > 0.6 fm~^, which 



is not consistent with microscopic calculations [25||. In this article, we shall propose new 
parameter-sets of the M3Y-type semi-realistic A^A^ interaction. As far as the energy of the 
symmetric nuclear matter is fixed, the symmetry energy at each density is well connected 
to the energy of the neutron matter. The new parameters are determined by fitting the 

(FP) or 



neutron-matter energy to microscopic result in Ref. |26 



[25] (APR). Moreover, 

we additionally take into consideration the binding energy of ^"'^Sn. As argued later, the 
symmetry energy at the saturation point tends to be fixed with good precision by fitting 
the parameters both to ^^''Sn and ^'^^Sn. The symmetry energy is thus constrained to good 
degree in the new parameter-sets. Corresponding to the microscopic results on the neutron- 
matter energy, we obtain two parameter-sets 'M3Y-P6' (fitted to FP) and 'M3Y-P7' (to 
APR). Although the parameters are determined from a limited number of data, they will 
be useful for investigating various aspects of nuclear properties, as will be illustrated by 
separation energies of the proton- and neutron-magic nuclei and by Z- or A^-dependence of 
shell structure. 
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II. M3Y-TYPE INTERACTION 



We take a non-relativistic isoscalar nuclear Hamiltonian of 

H^ = K + V^; K = J2^, V^ = J2^,,, (1) 

i i<j 

with i and j representing the indices of individual nucleons. We set M = (Mp + M„)/2 
throughout this paper, where Mp (M„) is the mass of a proton (a neutron) [23]. For the 
effective A^A^ interaction Vij, the following form is considered, 

(C) , (LS) , (TN) , (DD) 



n 

r = E i^'n'^^'PTE + e^°)PT0)/rHr.) L., ■ (S. + S,) , 

n 

= ihPs. ■ Ipir.r'^' + t^^TE ■ [Pir.rn Sir,) , (2) 



where Sj is the spin operator of the i-th nucleon, rij = — r^, r^j = |rjj|, pjj = (pj — Pj)/2, 
Lij = Tij X pij, and p(r) denotes the nucleon density. The tensor operator is defined by 
Sij = 4 [3(sj ■ fij)(sj ■ fij) — Sj ■ Sj] with f-jj = Vij/rij. The projection operators on the singlet- 
even (SE), triplet-even (TE), singlet-odd (SO) and triplet-odd (TO) two-particle states are 

1 - 1 + „ 1 + P, 1 - P, 



SE 2 2 ' 2 2 



_1-P. 1-P. _1 + P. 1 + P. 
Pso-^ ^To-^ (3) 

where P^- [Pr) expresses the spin (isospin) exchange operator. The Yukawa function /„(r) = 
e~'^'"^//i„r is assumed for all channels except v^^^\ The density-dependent contact term 
y(DD) jg g^fjfjgfj order to reproduce the saturation properties. Physically, t>(°°) may carry 
effects of the NNN interaction and of the density- dependence that is dropped in the original 
M3Y interaction. 

We start from the M3Y-Paris interaction 0, which will be denoted by M3Y-P0 in this 
article as in Ref. The range parameters /i„ of M3Y-P0 are maintained in any of v'''^\ 
v^^^^ and w*^"'"^''. As in M3Y-P0, the longest-range part in v^'^'^ is kept identical to the central 
channels of the one-pion exchange potential (OPEP), fopEP- Although the is splitting plays 
a significant role in the nuclear shell structure, the G-matrix is known to underestimate the 
is splitting. Even though effects beyond the MFA may cure this problem 28l| , we introduce 



an overall enhancement factor to v^^^^ in order to describe the shell structure within the 
MFA. Effects of the tensor force on the single-particle (s.p.) levels could be relevant to the 
new magic numbers in unstable nuclei js], We keep v^™'' without any modification from 
M3Y-P0. Because of this v^™'^ having realistic nature, the present M3Y-type interactions 
are useful to investigate the tensor-force effects within the MFA and RPA, as shown in 



Refs. |13|-|15| with the previous parameter-sets. The parameters in M3Y-P6 and P7 are 



tabulated in Table [U together with M3Y-P0. 
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TABLE I: Parameters of M3Y-type interactions. 



parameters 




M3Y-P0 


M3Y-P6 


M3Y-P7 




(fm) 


0.25 


0.25 


0.25 


^(SE) 


(MeV) 


11466. 


10766. 


10655. 


(TE) 
''1 


(MeV) 


13967. 


8474. 


9592. 


^(SO) 


(MeV) 


-1418. 


-728. 


11510. 


^(TO) 


(MeV) 


11345. 


12453. 


13507. 




(fm) 


0.40 


0.40 


0.40 


(SE) 
''2 


(MeV) 


-3556. 


-3520. 


-3556. 


(TE) 
2 


(MeV) 


-4594. 


-4594. 


-4594. 


(SO) 
2 


(MeV) 


950. 


1386. 


1283. 


(TO) 
2 


(MeV) 


-1900. 


-1588. 


-1812. 




(fm) 


1.414 


1.414 


1.414 


(SE) 
3 


(MeV) 


-10.463 


-10.463 


-10.463 


(TE) 


(MeV) 


-10.463 


-10.463 


-10.463 


(SO) 
3 


(MeV) 


31.389 


31.389 


31.389 


(TO) 
3 


(MeV) 


3.488 


3.488 


3.488 




(fm) 


0.25 


0.25 


0.25 


^(LSE) 


(MeV) 


-5101. 


-11222.2 


-11732.3 


^(LSO) 


(MeV) 


-1897. 


-4173.4 


-4363.1 




(fm) 


0.40 


0.40 


0.40 


(LSE) 
''2 


(MeV) 


-337. 


-741.4 


-775.1 


(LSO) 
''2 


(MeV) 


-632. 


-1390.4 


-1453.6 




(fm) 


0.40 


0.40 


0.40 


.(TNE) 
^1 


(MeV-fm-2) 


-1096. 


-1096. 


-1096. 


,(TNO) 
1 


(MeV-fm-2) 


244. 


244. 


244. 




(fm) 


0.70 


0.70 


0.70 


.(TNE) 
2 


(MeV-fm-2) 


-30.9 


-30.9 


-30.9 


,(TNO) 
2 


(MeV-fm-2) 


15.6 


15.6 


15.6 


a(SE) 







1 


1 


(SE) 


(MeV-fm^) 


0. 


384. 


830. 








1/3 


1/3 


(TE) 


(MeV-fm) 


0. 


1930. 


1478. 



III. PROPERTIES OF NUCLEAR MATTER 

We first apply the new semi-realistic interactions to the infinite nuclear matter in the 
Hartree-Fock (HF) approximation. Notice that only v^*"^ + in Eq. ([2]) contributes to 
the nuclear matter properties in the MFA. Energy of the nuclear matter is a function of the 
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following variables, 



P 



^Pra 

(TT 



Vt 



^ Ppt ~ Ppj + Pnt ~ Pnj 

P P 

Pp\ + Ppi - Pn\ - Pnl 



E 



P P 



^ ^ Ppt ~ Ppi - Pn^ + Pnj /^x 

where t = p,n and a =t, 4 are sometimes substituted by ±1 without confusion. As we 
restrict ourselves to the properties at zero temperature, the density depending on the spin 
and the isospin pra is related to the Fermi momentum kpra via 

Pr<T = '^^Fra ■ (5) 

Basic formulas to calculate the nuclear matter energy and its derivatives for given kpra were 
derived in Ref. Note that the superfluidity barely influences the nuclear matter energy, 
even if it takes place. 

We shall denote energy per nucleon {E/A) by S, where E is the expectation value of Hjsi 
for the nuclear matter. The spin-saturated symmetric matter is characterized by rjs = rjt = 
rjst = 0, for which we represent kpra simply by kp. Minimization of S{p), 



dS_ 

dp 



dkv 



0, (6) 



determines the saturation density po (equivalently, kpo) and the saturation energy Sq. The 
expression |o indicates evaluation at the saturation point. 

We depict £{p) for the spin-saturated symmetric nuclear matter in Fig. [H up to p ~ 5po. 
The results of the new semi-realistic {i.e. M3Y-P6 and P7) interactions are compared with 
those of the Skyrme SLy5 0, the Gogny DIS Q and DIM |3l| interactions. While all 



effective interactions give close S{p) at p < po, interaction-dependence is visible at p > 2po, 
though the SLy5 energy almost coincides with the M3Y-P6 one. 

Energy per nucleon in the spin-saturated neutron matter {i.e. rjt = —1, rjs = rjst = 0) is 



shown in Fig. [21 The FP [26j and APR [25[ results, to which M3Y-P6 and P7 are respectively 
fitted, are also presented. Having been fitted to a microscopic result [3^ as well, the energy 
with SLy5 is close to that with M3Y-P7 at any p. The stronger p dependence in M3Y-P6 
and P7 than in DIS and DIM originates from the choice a^^^^ = 1 in v^^^\ and enables us 
to reproduce the microscopic results. Since ti(^^) drives repulsion in the SE channel at high 
p but does not in the TO channel, the interactions having the form of Eq. ([2]) may give rise 
to the spin-polarized phase at high p in the pure neutron matter. However, the transition to 
the spin-polarized phase is delayed until p ~ 9po (20po) for M3Y-P7 (P6), almost irrelevant 
even to the neutron stars. 
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150- 




0.0 0.2 0.4 0.6 0.8 

p [fm"^] 



FIG. 1: (Color online) E = E/A vs. p in the symmetric nuclear matter, calculated with M3Y-P6 
(red dashed line), M3Y-P7 (red solid line), DIS (blue solid line), DIM (blue dashed line) and SLy5 
(green dot-dashed line). 




0.0 0.2 0.4 0.6 0.8 



p [fm"^] 

FIG. 2: (Color online) £ = E/A vs. p in the pure neutron matter. Circles and diamonds represent 
the FP and APR results, respectively. See Fig. [T] for the other conventions. 
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TABLE II: Nuclear matter properties at the saturation point. 







M3Y-P6 


M3Y-P7 


DIS 


DIM 


SLy5 




(fm-1) 


1.340 


1.340 


1.342 


1.346 


1.334 


So 


(MeV) 


-16.24 


-16.22 


-16.01 


-16.02 


-15.98 


/Co 


(MeV) 


239.7 


254.7 


202.9 


225.0 


229.9 


M^/M 




0.596 


0.589 


0.697 


0.746 


0.697 


am 


(MeV) 


32.14 


31.74 


31.12 


28.55 


32.03 


O-sO 


(MeV) 


26.47 


23.04 


26.18 


16.56 


37.47 


dstO 


(MeV) 


41.00 


43.30 


29.13 


28.71 


15.15 


Qo 


(MeV) 


-378.0 


-320.1 


-515.7 


-459.0 


-363.9 




(MeV) 


44.64 


51.53 


22.44 


24.83 


48.27 



The symmetry energy is defined by the second derivative of £ with respect to rjt for the 
spin-saturated matter, 

Id's 
at{p) = - 



(7) 



Here \p indicates rjt = rjs = rjst = but p left as a variable. The symmetry energy at the 
saturation point at{po) is denoted by a^o- Analogously, we define 



Id's 

as{p) = - 



Id's 
ast{p) = - 



2dr]i 



(8) 



and a^o 
by 



O'siPo), cisto = CLstipo)- The incompressibility at the saturation point is obtained 



/Cn 



B'S 
dkl 



9p' 



d'S 



(9) 



9p 

The effective mass (fc-mass) is defined by a derivative of the s.p. energy eihar). We denote 
the effective mass at the saturation point by Mq, which is given as 



de(k(7T) 



dk 



(10) 



These characteristic quantities calculated from the new semi-realistic interactions are tab- 
ulated in Tables [Tll Those from DIS, DIM and SLy5 are also displayed for comparison. 
Also compare them with the empirical values fcpo ~ 1-33 — 1.34 fm~^, Sq ~ —16 MeV, 



33| and ato ~ 30 MeV l34|- It is remarked that the fit both to ""Sn and 



/C ~ 240 MeV 

"^^^Sn (see Table llVjl well constrains ato in the M3Y-type interactions with good precision, 
to ato ~ 32 MeV. This ato value is in harmony with ato of SLy5, though not so wit h am 



of DIM and of a recently proposed Skyrme-type energy density functional UNDEFl [35 



The effective mass (Mq ^ 0.6 M) of the current M3Y-type interactions is not much change- 
able, which does not contradict to a microscopic result |36| but is lower than the value that 
reproduces collective excitations in the RPA {e.g. the DIM value). 
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As mentioned in Introduction, p-dependence of the symmetry energy attracts interest. 
The first derivative of at{p) at po is under debate, which is parametrized as 



d 

Ao = 3 -rat[p) 
dp 



2 ^ dkp drj- 



3 d'^S 
2^ dpd-q 



(11) 



The characteristic coefficient £to; along with the third derivative of £ with respect to p that 
is denoted by Qo; 



Qo 



27p^ 



(12) 



_ 

are also presented in Table HIl It is noteworthy that M3Y-P6 and P7 have higher Cto than 
DIS and DIM, in contrast to the previous parameter-set M3Y-P5 ^ that has comparable 
Cto to DIS and DIM. The higher values seem favorable for describing the low-lying El 
strengths [23 . 

The symmetry energy at{p) in a wider region of p is depicted in Fig. [3l along with 
as{p) and ast{p)- If any of at(p), a,s(p) or ast{p) is negative, the spin-saturated symmetric 
nuclear matter becomes unstable, undergoing phase transition. With DIS the symmetric 
matter is unstable beyond p ~ 3.4po, as inferred from Fig. [2] and manifested in Fig. [3]-a). 
Moreover, Fig. [3]-c) implies that the magnetized phase emerges at moderately high p, when 
we employ SLy5 or DIM. The transition takes place at p ~ 2.1po (3.0po) in the SLy5 (DIM) 
result. Similar instability occurs at p ~ (1.2 — 3)po for most available parameter-sets of the 
Skyrme interaction including the tensor channels, even if deformation of the Fermi sphere 



is ignored [37j. On the contrary, in this density region the isotropic nuclear matter is stable 
against spin or isospin asymmetry under M3Y-P6 and P7, though M3Y-P6 gives decreasing 
at{p) in p > 2.2po, which eventually becomes negative in p > 5.8po. 

The Landau-Migdal (LM) parameters have been used to argue global characters of ex- 
citation modes of nuclei. Employing the analytic formulas given in Ref. [3], we evaluate 
the LM parameters at the saturation point for the new semi-realistic interactions, as shown 
in Table lllli See Ref. 0] for definition of the LM parameters. Several LM parameters are 
related to the characteristic coefficients in Table [III as follows. 



Ml 
M 



1 + 3A 



/Co 



3P 
■ Mo* 

^FO 



I + /0) 



1 + 



6M0* 



;i + /o) 



(13) 



It has been known that go is small while is relatively large (~ 1) [38|. M3Y-P6 and 
P7 hold reasonable characters on the spin and isospin channels as the previous parameters, 
owing significantly to f qpep 



IV. APPLICATIONS TO FINITE NUCLEI 

We next apply the new semi-realistic interactions to finite nuclei within the MFA. The 
Hamiltonian H = + Vc — -ffc.m. is used, with given in Eq. ([T]). Vc and i^c.m. 
represent the Coulomb interaction and the center-of-mass (cm.) Hamiltonian. We have 
made no additional approximation on H, by handling the exchange term of Vc and the 
two-body terms of i^c.m. explicitly. Effects of Vc on the proton pairing, which have recently 
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TABLE III: Landau-Migdal parameters at the saturation point. 





M3Y-P6 


M3Y-P7 


DIS 


DIM 


SLy5 


/o 


-0.360 


-0.329 


-0.369 


-0.255 


-0.276 


fl 


-1.211 


-1.233 


-0.909 


-0.762 


-0.909 




-0.394 


-0.381 


-0.558 


-0.302 


0.0 


h 


-0.183 


-0.177 


-0.157 


-0.058 


0.0 


fo 


0.544 


0.506 


0.743 


0.701 


0.815 




0.511 


0.571 


0.470 


0.378 


-0.387 




0.225 


0.234 


0.342 


0.632 


0.0 


fs 


0.090 


0.091 


0.100 


0.137 


0.0 


90 


0.272 


0.093 


0.466 


-0.013 


1.123 


91 


0.231 


0.337 


-0.184 


-0.380 


0.253 


92 


0.163 


0.179 


0.245 


0.483 


0.0 


93 


0.077 


0.079 


0.091 


0.114 


0.0 


9'o 


0.970 


1.055 


0.631 


0.711 


-0.141 


9'i 


0.157 


0.069 


0.610 


0.652 


1.043 


9'2 


0.053 


0.044 


-0.038 


-0.243 


0.0 


9'3 


0.004 


0.004 


-0.036 


-0.064 


0.0 



been recognized to be sizable [39[, are explicitly included in the Hartree-Fock-Bogolyubov 
(HFB) calculations as well. 

The algorithm based on the Gaussian expansion method (GEM) j^, loj is applied for all 
the numerical calculations of finite nuclei in this article. In this method we employ the s.p. 
basis-functions of 



uij \ 



(14) 



apart from the isospin index, where Y^^\r) is the spherical harmonics and Xa- the spin wave 
function. For the range parameter u, which is generally a complex number [u = Uj- + iui), 
we adopt the following values lo| : 



-2n 



±- 



TT 



[n 
in 



0,1,--- ,5) 

0,1,2) 



(15) 



with z/q = (2.40 fm)~^ and h = 1.25, resulting in 12 bases for each In the HFB 

calculations the s.p. space is truncated as ^ < 8 (for the Z = 82 and = 126 nuclei) or 
i < 7 (for the lighter ones). As shown in Ref. 10|, the above set of the GEM bases can 
cover wide range of nuclear mass with good precision. 

In the following we shall show results of the spherical HF and HFB calculations using 
the new M3Y-P6 and P7 interactions, in comparison with those using the Gogny DIS 
and DIM interactions. Although there have been more advanced calculations using the 
Gogny interactions, e.g. the calculations based on the generator-coordinate method jZoj, 
comparison is made only at the MFA level in this paper, leaving extensive applications of 
the M3Y-type interactions as a future work. 
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0.0 0.1 0.2 0.3 0.4 0.5 0.6 
p [fm"'] 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 

p [fm"'] 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 



p [fm"'] 

FIG. 3: (Color online) a) at{p), b) as{p) and c) ast{p) in the symmetric nuclear matter. See Fig. [1] 
for the conventions. 
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TABLE IV: Binding energies {—E) and rms matter radii (a/ {r'^)) of several doubly magic nuclei. 
Experimental data are taken from Refs. [4l|444i|. 









Exp. 


M3Y-P6 


M3Y-P7 


DIS 


DIM 


160 


-E 


(MeV) 


127.6 


126.3 


125.9 


129.5 


128.2 






(fm) 


2.61 


2.59 


2.57 


2.61 


2.57 


240 


-E 


(MeV) 


168.5 


166.2 


167.4 


168.6 


167.3 




V \ ' 


(fm) 


3.19 


3.05 


3.03 


3.01 


2.98 


40Ca 




(MeV) 


342.1 


335.9 


334.3 


344.6 


342.2 




\/{r'^) 

V \ ' 


(fm) 


3.47 


3.37 


3.35 


3.37 


3.33 


48Ca 




(MeV) 


416.0 


413.8 


414.9 


416.8 


414.6 




V \ / 


(fm) 


3.57 


3.51 


3.49 


3.51 


3.48 


902r 


-E 


(MeV) 


783.9 


781.1 


780.8 


785.9 


782.1 




Vir') 


(fm) 


4.32 


4.23 


4.22 


4.24 


4.20 


lOOSn 


-E 


(MeV) 


824.8 


822.5 


822.8 


831.6 


824.9 






(fm) 




4.36 


4.34 


4.36 


4.32 


i32Sn 


-E 


(MeV) 


1102.9 


1097.8 


1100.8 


1104.1 


1104.5 






(fm) 




4.78 


4.77 


4.77 


4.72 


208pb 


-E 


(MeV) 


1636.4 


1634.5 


1635.5 


1639.0 


1638.9 






(fm) 


5.49 


5.53 


5.51 


5.51 


5.47 



A. Doubly magic nuclei 

The spherical HF approach is rationally expected to be a good approximation for the 
ground states of the doubly magic nuclei. We present the binding energies and the rms matter 
radii of several doubly magic nuclei in Table IIVI The spherical HF results using the new 
semi-realistic interactions are compared with those using DIS and DIM, as well as with the 
experimental data. Influence of the cm. motion on the matter radii is subtracted in a similar 
manner to the cm. energies [4|. The binding energies of these nuclei by the M3Y-P6 and 
P7 interactions are in agreement with the measured values within 5 MeV accuracy, except 
^°Ca. This accuracy is comparable to those of DIS and DIM. In ^°Ca, influence of octupole 
correlations might be strong, as suggested by the low 3^ energy in measurements jisi and 
mentioned in Ref. [5|. For this reason we have not taken this discrepancy seriously at the 
present stage, while future study is needed on this problem. The rms matter radii of these 
nuclei calculated from M3Y-P6 and P7 are also in fair agreement with the data. We point 
out that DIS has not predicted accurate energy of ^"^Sn, and that DIM systematically gives 
smaller radii than the measured ones. 

The non-central channels of the effective interaction, v^^^^ and f are responsible for 
the is splitting of the s.p. levels and its nucleus-dependence. We display the s.p. levels of 
208p|j calculated in the HF approximation, comparing to the observed levels in Fig. HI The 
experimental s.p. levels are taken from the lowest states having specific spin-parity in the 
A = 207 or 209 nuclei. Because of the fragmentation via the coupling to the many-particle- 
many-hole configurations, these observed states do not straightforwardly correspond to the 
s.p. levels in the MFA. In M3Y-P6 and P7 the non-central channels are not changed from 
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FIG. 4: Single-particle energies for ^'^^Pb. Experimental values are extracted from Refs. [4ll. l45l|. 



TABLE V: Difference between proton and neutron rms radii y^(r^ - in ^°^Pb (fm) by 

the HF calculations with several interactions. 



M3Y-P6 


M3Y-P7 


DIS 


DIM 


0.158 


0.161 


0.136 


0.112 



M3Y-P0 except the overall enhancement factor to v^^'^\ This factor is determined so that 
the level ordering should not differ seriously from the observed one around ^"^^Pb. While 
appropriateness of the enhancement factor to should further be investigated in future 
studies, it is a simple and useful cure to the Is splitting. 

It has been recognized that the neutron-skin thickness in ^osp]^ jg connected to the p 
dependence of the symmetry energy 4^ , particularly the Cto parameter. Table |V] presents 
difference between proton and neutron rms radii \J {r'^)n — \/Jt^)p in ^°^Pb. The values of 
M3Y-P6 and P7 are in good agreement with the recent experimental value drawn from the 
El strengths [47| . 



B. Proton- or neutron-magic nuclei 



The spherical HFB approach provides us with a reasonable approximation for the nuclei 
in which Z or is a magic number. The odd-Z or nuclei can be handled in the equal- 
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FIG. 5: (Color online) Neutron separation energies for a) Z = 8, h) Z = 20, c) Z = 28, d) Z = 50 
and e) Z = 82 nuclei, calculated with M3Y-P7 (red circles) and DIM (blue open triangles). Lines 
are drawn to guide eyes. Experimental values are taken from Ref. [41] and presented by the crosses. 



filling approximation [48|, |49[. In fixing the new parameter-sets, we have taken into account 
the pairing properties by fitting to the data on the even-odd mass differences in the Z = 50, 
~ 70 and the = 82, Z ~ 60 nuclei 0. 

In Fig. |5] (Fig. [6]), the neutron (proton) separation energies S'„ (Sp) are plotted for the 
Z = magic (A^ = magic) nuclei. The Sn and Sp values of M3Y-P6, which are always close to 
those of M3Y-P7, are not presented. Similarly, being close to the DIM ones, the DIS results 
are not displayed. Notice that the even-odd mass difference is proportional to the difference 
of the separation energies between the adjacent nuclei, while the two-neutron (two-proton) 
separation energy is the sum of SnS (5'p's) of the two neighboring nuclei. Although there are 
certain discrepancies if we look into their details, M3Y-P6 and P7 give separation energies 
in agreement with the measured ones with the accuracy similar to DIS and DIM. 

It has been suggested that isotope shifts of the Pb nuclei may be relevant to the is 
potential [sT], which is primarily determined by v^^^K In 6^y {r'^)p{'^Ph) = a/ (r^)p('^Pb) — 
V (^^)p(^°'^Pt)), a kink has been observed at ^°*Pb. The zero-range LS force contained in the 
Gogny as well as the original Skyrme interactions operates only on the T = 1 two-nucleon 
states. It was argued that this isospin character of the LS force could be insufficient to 
describe the rapid rise of 6^y {r'^)p in A^ > 126. Having v'-^^'^ with finite ranges that acts also 
on the T = channel, it may be interesting whether the M3Y-type interactions reproduce 
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the kink of 5^J {r^)p- 

In Fig. [71 we depict ^ a/ {r'^)p{'^Ph) obtained by the spherical HFB calculations with 



53 . The results of 



M3Y-P7, DIS and DIM, in comparison with experimental data [52|, 
M3Y-P6 are almost indistinguishable from those of DIM in iV < 126 and from those of 
M3Y-P7 in > 126. The DIS interaction does not give clear bending at ^°^Pb. On 
the contrary, although its LS force holds the zero-range form, DIM provides a visible kink 
at 208pb. M3Y-P7, in which v^^^^ has finite ranges, further improves 6^y {r'^)p. We have 
confirmed that this tendency well correlates to the occupation probability on nOiii/2, as 
pointed out in Ref. js^. It seems reasonable to consider that the s.p. energy difference 
e(n0iii/2) — e{nlgg/2) is responsible for the interaction-dependence of S^y {r^)p in > 126. 
However, whereas isospin character of t>('"^) plays a certain role in the s.p. energy difference, 
it is not yet obvious whether and how interplay of other channels, e.g. v^™^ and the pairing, 
contributes to 6 

In earlier studies 
Ref. M 




the MFA results were compared with the data given in, e.g., 
The kink at ^osp]^ becomes stronger in the new data [s^ because of larger values 



of 5^^/ {r'^)p in > 126, which has not been reproduced within the MFA, to the author's 
best knowledge. 
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FIG. 7: (Color online) Isotope shifts of the Pb nuclei 5 ^JJr^Yp{^Ph) , obtained from the HFB 
calculations with M3Y-P7 (red solid line), DIM (blue dashed line) and DIS (brown dot-dashed 



line). Experimental data are taken from Refs. 
crosses, respectively. 



and 



53], which are represented by squares and 



C. Shell structure and tensor force 



It has been established that the shell structure depends on the effective A^A^ interaction of 
the MFA. In particular, the effects of the tensor force are significant, as clarified in the proton- 
or neutron- magic nuclei [sl, 17|. In some nuclei the spin-isospin channels of the central force 
could also have appreciable effects jl, [l6|. Having realistic v^™'' and fopEP; M3Y-P6 and P7 
are suitable to investigating those effects as the previous interaction M3Y-P5 [1, [l8|. This is 
a clear advantage of these semi-realistic interactions over the phenomenological interactions 
such as DIS and DIM. Although the results of M3Y-P6 and P7 on the shell structure are 
not essentially different from those of the previous parameter-set shown in Refs. [sl, [13, 18|. 
we briefly mention several notable points. 

The iV-dependence of the differences of the observed s.p. energies e{p0hii/2) — £{pldr,/2) 
and e{p0gj/2) —e{pld5/2) in the Sn isotopes are simultaneously reproduced, with being 
crucial to the former (see Fig. 14 and related arguments in Ref. |5|). In the = 32 isotones 
y(TN) v^p-^p give rise to significant Z-dependence of ^(^0/5/2) — £(^1^3/2) (see Fig. 13 
of Ref. [sl and Fig. 3 of Ref. [l3]), possibly accounting for the new magic number = 32 
around ^^Ca 55|, |56j. However, despite the presence of the realistic tensor force, the semi- 
realistic interactions do not predict closure of A^ = 34 at ^^Ca (see Fig. 3 of Ref. [isj), in 



contrast to the shell model prediction in Ref. |57|. While semi- magic nature of A^ = 40 
is indicated at ^^Ni with the semi-realistic interactions as with the Gogny interactions (see 
Fig. 4 of Ref. 18|), which seems consistent with experiments [58|, it is likely for the A^ = 40 
magic nature to be broken at ^°Ca because of v^™^ (see Fig. 3 of Ref. [l8| and Fig. 4 of 
Ref. [1^). It is of interest to investigate magicity of A^ = 40 toward ^°Ca experimentally, 
which could further clarify role of f The doubly magic nature of ^^Ni is predicted to 
hold even with v^™\ giving ^^(2+) = 3.0 - 3.5 MeV and B{E2; 2+ 0^) ^ 85 eW in the 
HF-fRPA. Experimental data on 2f of ''^Ni are awaited. 

Location of the neutron drip line for the Ca and Ni nuclei could be investigated by the 



new experimental facilities [59j and has been argued in Ref. [l4j. We tabulate the neutron 
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TABLE VI: Neutron numbers of the heaviest bound Ca and Ni nuclei predicted by the spherical 
HFB calculations with several interactions. 



Isotope 


M3Y-P6 


M3Y-P7 


DIS 


DIM 


Ca 


50 


50 


44 


46 


Ni 


60 


62 


58 


58 



drip line predicted by the spherical HFB calculations with the M3Y-type and the Gogny 
interactions, in Table IVII 

V. SUMMARY AND OUTLOOK 

We have developed new parameter-sets of the semi-realistic effective interactions to de- 
scribe low energy phenomena of nuclei. They are obtained by phenomenologically modifying 
several parameters in the M3Y-Paris interaction, while the tensor force and the OPEP part 
in the central force are not changed, as before. Unlike the previous parameters, the new 
sets M3Y-P6 and P7 are adjusted also to the microscopic (FP and APR) results of the 
neutron-matter energies and to the binding energy of ^°'^Sn. We therefore attain improve- 
ment on the symmetry energy, up to its density-dependence. In contrast to instability of 
the spin-saturated symmetric nuclear matter in the SLy5, DIS and DIM results, neither of 
M3Y-P6 nor P7 predicts such phase transition in the density range of p < 6po- 

The new parameter-sets M3Y-P6 and P7 have been applied to the doubly magic nuclei 
in the spherical HF calculations, and to the proton- or neutron-magic nuclei in the spherical 
HFB calculations. Fair agreement with experimental data has been demonstrated for the 
binding energies of the doubly magic nuclei and for the nucleon separation energies of the 
proton- or neutron-magic nuclei. Owing to the realistic tensor force and the OPEP central 
force, the Z- or iV-dependence of the shell structure is well described with M3Y-P6 and P7, 
as with the previous set M3Y-P5. The isotope shifts of the Pb nuclei have also been argued. 

Future study includes application of the semi-realistic interactions to the excitations in 
the RPA, as well as to deformed nuclei. Moreover, extensive applications to nuclear reactions 
and to the neutron stars may be within the reach, which give further test of the effective 
interactions and a step toward unified description of nuclear structure, reactions and neutron 
stars. 
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